The Random Fuse Network as a Dipolar Magnet 
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We introduce an approximate mapping between the random fuse network (RFN) and a random 
field dipolar Ising model (RFDIM) . The state of the network damage is associated with a metastable 
spin configuration. A mean-field treatment, numerical solutions, and heuristic arguments support 
the broad validity of the approximation and yield a generic phase diagram. At low disorder, the 
growth of a single unstable 'crack' leads to an abrupt global failure. Beyond a critical disorder, the 
conducting network sustains significant damage before the coalescence of cracks results in global 
failure. 



Heterogeneities are present in real materials, in a wide range of magnitudes and scales, and generate fluctuations 
in local toughness. The fracture mechanics then results from an interplay between disorder and stress fluctuations 
jl],^| . In attempting to elucidate the role of the disorder, a number of classifications and characterizations of rupture 
have been proposed |3|-|l0| . In particular, the behavior of the damage as the applied force approaches its critical value 
(global rupture) — partial versus no damage, macroscopic abrupt versus continuously increasing damage — allowed to 
distinguish various types of fracture. Most studies were carried out either numerically or in the framework of extreme 
statistics (see Refs. JUHI and references therein), while few satisfactory analytical approaches exist. In parallel, a 
number of authors jll],[12| introduced the so-called random fuse network (RFN), less involved than fracture models 
because of its scalar nature but still capturing the main aspects of rupture mechanics |l3| . Thus the RFN often is 
seen as a toy model for the investigation of mechanical rupture. The mechanical bonds are replaced by fuses burning 
if the current (or equivalently the electrical field) exceeds a random threshold. As the total applied current flowing 
through the network is increased, more and more fuses burn, thereby decreasing the effective conductivity of the 
network until a global failure at which no current can go through anymore. The similarity of the RFN with a driven 
random field Ising model at zero temperature ]14| ] has been recently noted || , but no systematic correspondence has 
been established. 

In this letter, we present a mapping, valid at least at the initial stages of failure, between the RFN and a driven 
random field Ising model with dipolar coupling. This mapping associates the configuration of the burnt fuses at a given 
value of the applied electric field with a zero temperature metastable state of the random field dipolar Ising model 
(RFDIM). The mapping between fracture and a classical model of statistical mechanics, is supported by mean-field 
and heuristic arguments and by a preliminary numerical study, from which a generic phase diagram with a disorder 
induced transition emerges. 

The RFN consists of a d-dimensional lattice of fuses, each with a conductivity ui (x) = oq initially, where x = am, 
m G Z d and a is the lattice spacing (we will assume the continuum limit a — ► 0), and i = 1, . . . , d corresponds to the 
orientation of the fuse. If the electric field along i at x, 2£j(x), exceeds a random threshold £j(x), the fuse burns and 
<7j(x) = 0. We consider a scenario in which the electric field in the direction i = 1 is slowly increased from zero and 
Ej = for % ^ 1 at all times. 

Assuming locality and analyticity, we can write a general equation for the time evolution of the conductivity field 

as 



where Aj = £i(x) 2 — £^(x, t) 2 and the Heaviside function ensures that the conductivity of the fuse remains Co as 
long as the electric field it sustains does not exceed its threshold. The simplest choice for the function /j within the 
constraints of positivity ((7j(x, t) > 0) and monotonicity (9t(jj(x, t) < 0), /, = — constant x <Tj(x, t), yields 



after absorption of the constant in a redefinition of r\. Equation (Q) incorporates two time scales: the relaxation 
time r\ and the characteristic time r = E\{dE\ j d£)~ x . We focus on the limit r\jr — » 0, in which the conductivity 
relaxes instantaneously from <tq to as the local field crosses the local threshold. To capture the sequential dynamics 
associated with the r\ ^ retardation effect, we discretize |JL5) Eq. (0) in time as 



77cVi(x,t) = f i ({a j {x,t)},{d j a k (x,t)},...)0(-A i (x.,t)), 



(1) 



CTi(x,i)0(-Ai(x,t)) 



(2) 



Oi (x, t + St) 



£7 i (x,t)0(A i (x,t)). 



(3) 
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Formally, this corresponds to setting rj = St and discretizing Eq. (0) in the usual fashion. As expected, Eq. (||) 
prescribes an irreversible dynamics such that <7i(x) = at time t implies <Ti(x) = for all later times. 

For a given configuration of the conductivities, Kirchoff's laws may be recast [|l6| into the integral equation 

£i(x) = E 5 a + f tfy^Gyfx - y)Saj(y)Ej(y), (4) 
J j 

where Eq is the magnitude of the external field applied along the first axis, (5<7j(x) = <7j(x) — <tq, G,j (x) = —didjg(x) 
the dipolar tensor [17a], and the Green's function g is defined by croV 2 g(x) = — <5(x). To lowest non-trivial order, 
Kirchoff's laws become 



Ei(x) = E Q 



J d d yG il (ic-y)6a 1 (y) 



(5) 



As we shall see, this simple approximation captures much of the physics of network damage. First, due to the presence 
of the dipolar kernel Cry, Eq. (@) reproduces the expected shielding of the current 'upstream' and 'downstream' to a 
burnt region, as well as the lateral amplification Second, the fcth term of the expansion contains a product 

Ili=i fc^ cr i(y«) which vanishes with high probability, especially at the early stages of the damage. In addition, this 
approximation satisfies the necessary property that the field at the tip of a crack increases with crack growth. For one 
isolated burnt fuse, the field on the next parallel link in the direction perpendicular to the applied field is enhanced 
by a factor a > 1. The approximation (||) gives a = 1 — Gn(0, 1) ~ 2.14 on the 2d square lattice (and a = 2 in the 
continuum limit [17b]); the exact result on the square lattice is a — 4/tt. When n such parallel fuses are broken, the 
field at the tip of the crack is known to scale as n 1 / 2 ^ -1 ) while in the approximation of Eq. (||) 

E tip = E 1 (0,p = n + l) = E o (l+ J2 X /P d ) (6) 

\<p<n+l 

(where p 2 = x\ + . . . + x\ is the transverse square distance) converges to a finite value as n increases, corresponding 
formally to the d = oo case. For uniform disorder distributed according to p{£) — 6(£ 2 — + w — £ 2 )/w 

with w < 2 or a binary disorder E t i p can be larger than the largest threshold and we expect the same results as if 
Etip ~ n 1 / 2 ^ -1 ). For a disorder distribution allowing with a high probability very large values of thresholds (such 
as distributions with fat tails), the scaling of Ea p with the size of the crack becomes relevant and we cannot use the 
approximation (^). Here, we present results for the uniform disorder while similar results are obtained for binary 
disorder. 

Squaring Eq. (^) we obtain, to lowest nontrivial order, the distance A, from the local threshold in terms of the 
damage 5a i, as 

&x{^t)~£ l (xf-El-2Elj d d yG 11 ^-y)8a 1 {y,t) (7) 



A j (x,i)^^(x) 2 - J B 2 



j d d y G jl (x-y)Sa 1 (y,t) 



(8) 



(with j y£ 1) which, along with the update rule of Eq. (||) completely determines the state of the network under an 
applied field E . 

In terms of the 'Ising spins' defined by s(x, t) — 2Sa% (x, t)/cr + 1, Eq. (||) simply stipulates that each spin s(x, t) 
will lie along a local field Ai(x, t). Eq. (Q) shows that this local field Ai results from a uniform applied component, a 
random component, and a spin-spin interaction mediated by the component G\\ of the dipolar kernel. Therefore, as 
the applied field is switched from to Eq, spins flip in avalanches until they reach the state s(x) = sgn[Ai(x, t — oo)] 
which minimizes the Hamiltonian 

U = -\J d d xd d ys(x)J(x - y)s(y) - J d d xh(x) S (x), (9) 

where J(x) = -E 2 a G n (x), h = H + £ 1 {x) 2 , and H = -E 2 [l - er Q G n (0)] = —Eq(1 + l/d). As the driving field \H\ 
is increased, more and more spins flip downward, and the state of the RFN follows the metastable (nonequilibrium) 
spin state of Ti, connected by its history to the initial condition s(x) = +1 for all x. 

While similarities between the RFN and the random field Ising model (RFIM) have been pointed out in the past, 
here a precise mapping is elucidated. We stress that, at odds with previous studies on the RFIM [pfj in which the 
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spin-spin coupling is taken to be short-ranged, constant, and isotropic, the coupling J(x) is rather peculiar as its 
magnitude is proportional to the applied field H and its modulation is dipolar. Even in the absence of disorder, the 
dipolar interaction in Ising spin models leads to a rich behavior ^2),^) n °t yet fully understood. In particular, a 
striped phase was identified |24|,^5| with alternating up and down spins. A similar scenario occurs in the RFN at very 
low (narrow) disorder. All the fuses remain intact until Eq reaches the lowest threshold and the corresponding fuse 
burns. The dipolar coupling then relieves its longitudinal neighbors while further stressing its transverse neighbors, 
who burn in turn, and a stripe of burnt fuses develops, eventually spanning the whole system (the network then stops 
conducting and the formation of additional stripes is consequently prohibited). At high (wide) disorder, a resilient 
fuse in the passage of a burnt stripe will stop the growth of the crack, and new cracks will nucleate elsewhere. This 
suggest a transition, upon increase of the disorder J7|,p[p^], from an abrupt regime characterized by the growth of 
a single macroscopic unstable crack to a continuous regime in which the system is significantly damaged before its 
global failure. 

For a qualitative understanding of the role of disorder, we first revert to a simpler version of the problem in which 
the dipolar kernel is replaced by a uniform infinite-range coupling J(x) = EqJ/N, where J > measures the strength 
of the interaction and N is the total number of spins (equivalently, fuses). This 'democratic' model is mean-field like 
in that the evolution of the local field (Eq. (^)) acting on a spin s(x) may be described in terms of a single degree of 
freedom, the fraction n_ of down spins (burnt fuses), 

Ai(x) = £i(x) 2 - El - 2JEln-, (10) 

and depends on the position only through the local random threshold. Upon increase of the applied field a number of 
spins flip, thus increasing n_ which in turn causes more flips, and so on until this avalanching terminates at n_(i?o), 
the metastable fraction of downward spins at an applied field Eq. We follow the graphical scheme of Ref. (27) to obtain 
the resulting phase diagram. For the uniform disorder of width w, clearly rt_ = as long as Eq < £^ n , the lowest 
threshold. If w < 2J£^, rt_ then jumps to 1 abruptly; if w > 2J£f n , on the other hand, n_ = (Eq — £f n )/(w — 2JEq) 
increases continuously with E up to n_ = 1 at Eq = (£^ + w)/(l + 2J). As expected, a critical width wq = 2J£f n 
of the disorder separates two regimes: at low disorder, a huge avalanche is leading to global failures, while for large 
disorder, tiny avalanches can be stable, until cracks coalesce. We note the broad validity of this picture, at least in 
the mean field. In contrast with the 'democratic fiber bundle model' pq , ^9[ for which it was shown |^7j that abrupt 
failures are but an artifact of a large discontinuity in the threshold distribution, here the phase diagram extends to 
distributions that are not continuous, e.g., to any uniform distribution on a support [£^,£^ + w}. These mean-field 
results are also recovered from the homogeneous saddle point of the partition function 

z= Yl e ~ m ( n ) 

Mx)=±i} 

This saddle point corresponds to the stability condition for a microcrack: if the electric field at the tip of the crack 
is larger than the largest threshold, the crack is unstable. This implies the existence of the critical width wq found 
above. It is quite remarkable that this result — which was already found by probabilistic methods jjj — appears here 
naturally as the uniform saddle point or mean field. 

Because of the form of the coupling in this simplified democratic model, we can carry out a similar investigation 
without recourse to the linear approximation of Eq. (^), but using Eq. (^) directly. Assuming J < 1/2 (without which 
the problem is ill-defined), we find not only a similar phase diagram but also, surprisingly, the same value for the 
critical disorder wq — 2J£^. At least in this mean field version of the original problem, the linear approximation 
is legitimate. This also suggests that when Eq. (||) is valid, the early stages of the damage dominate the whole 
breakdown process. 

Evidently, a mean field phase diagram may be significantly modified by fluctuations governed by a dipolar coupling. 
On the one hand the dipolar kernel decays as l/r d , on the other hand its angular dependence invests the model with 
a so-called nonmonotonicity: burning a fuse relieves the current on a fraction of neighboring fuses. Monotonicity 
is a simplifying feature of many driven systems, and its lack clearly introduces complications and a possibly richer 



phenomenology 26 1. Furthermore, while a mean field treatment predicts the evolution of damage, it fails to capture 
the random fluctuations of the breakdown field E^ as well as finite-size effects. These two features are however present 
in the approximation (|j). In particular, it is easily shown that the maximum of the electric field scales as InL, as 
a consequence of the long range of the dipolar tensor. Investigation of finite size effects in the framework of this 
approximation present a complementary route to most of the approaches to finite size effects, which rely on extreme 
statistics 130] and neglect correlations. 
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For a better grasp of the full problem we obtained numerical solutions of the RFN (using Eqs. (||4|)) and the 
RFDIM (using the approximation Eq. ([?])). We found (Fig. 1) that the RFDIM reproduces quite well the behavior 
of the RFN. In both cases, we observe a critical value wq of the disorder beyond which there is a large window of 
damage that preceeds the global failure. 

Before concluding, we propose a first attempt to describe the stability of cracks in the spirit of the Imry-Ma 
argument |3l],|22|]. In three dimensions, the creation of a spherical crack of radius R modifies the total energy (cf. 
Eq. ®) by 

AH w 2hR 3 - wR 3/2 + gR 3 (12) 

where the last term comes from the dipolar interaction (g > 0) and merely renormalizes the average of the random 
field [^2|. If the applied field Eq is small, the average h is positive and the formation of a large crack is prohibited. 
In this case, the typical size of a crack will be i?o ~ [w/ (2ft + g)] 2 ^ 3 - On the other hand, if 2h + g < it is favorable 
for a crack to grow indefinitely. A more precise argument based on an oblate ('penny-shaped') crack perpendicular to 
the applied field leaves the above conclusions unchanged. As w is increased, the typical size Ro of damaged regions 
becomes larger, up to a critical value wq at which the latter percolate. This occurs when Ro becomes comparable to 
the typical distance between nucleation sites, a[w/{H — E 2 n )] 1 ^ 3 for an uniform disorder. As long as w < wq, breakdown 
results from a single crack that spans the system when 2h + g becomes negative upon increase of the applied field H . 
When w > wq, localized cracks grow and percolate before 2h + g changes sign, which allows significant damage while 
the network is still conducting. This heuristic argument also substantifies the broad validity of the approximation of 
Eq. (||) a priori only valid at the initial stages of the damage: the breakdown is controlled either by an instability or 
by the coalescence of (small) cracks. 

In summary, we have presented an approximate mapping of the RFN to the RFDIM, and argued for its broad 
validity. The state of the RFN — in particular the extent of damage (number of burn fuses) — is mapped into a 
metastable spin state. From a mean-field investigation, numerical solutions, and heuristic arguments, a generic 
picture emerges, characterized by a critical value of the disorder. At low disorder, a macrosocopic crack grows until it 
spans the system. At large disorder, threshold fluctuations stabilize (micro)cracks and a significant precursor damage 
develops before the latter coalesce into a percolating structure. Our investigation were carried out for a bounded 
disorder. By contrast, a slowly decaying threshold distribution may modify the character of the transition or suppress 
it altogether, an interesting question for future study. 
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FIG. 1. Effective conductivity E versus Eo for uniform dis- 
order for two values of the disorder width w = 0.5 < wo and 
w — 1.8 > wo (for a system size 50 x 50 and averaged over 
100 configurations), (a) Calculation for the RFN by solving 
Eqs. (3,4). (b) Calculation for the RFDIM obtained from the 
numerical solution of Eqs. (6a, b). In both cases (a) and (b), 
for w < wo, the breakdown is abrupt and there is no fluctua- 
tion in the breakdown field (for (a), it is not fully abrupt, due 
to finite-size effects). For w > wo, the breakdown field fluc- 
tuates and there is some damage before complete breakdown 
of the system. 
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